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Abstract. We present a near-IR spectroscopic analysis of the intriguing DO white dwarf RE J0503-285. The IR spectrum fails 
to reveal evidence of the presence of a spatially unresolved, cool, late-type companion. Hence we have placed an approximate 
limit on the spectral-type and mass of a putative companion (later than M8, M< 0.085M o ). This result rules out ongoing 
interaction between the white dwarf and a close companion with M> 0.085M o as responsible for the discrepancies between 
the observed photospheric abundances and model predictions. As the possibility remains that there is a cooler companion lying 
beyond the detection threshold of this study we use our modelling to estimate the improvement in sensitivity offered by a Spitzer 
observation. 
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1. Introduction 

In accord with a theoretical prediction made several decades 
ago (Schatzman 1958 1, the atmospheres of white dwarfs are 
observed to be dominated either by hydrogen (DAs) or by he- 
lium (DO/DBs). Unimpeded, their high surface gravities lead 
to the settling out of heavier elements on timescales of mere 
days to months (e.g. Dupuis et al. 1993 1. In the standard theory 
of single star evolution the distinction between the two com- 
positions is made at birth. If a star leaves the AGB during a 
period of quiescent hydrogen shell burning it evolves onto the 
hydrogen dominated cooling channel. Alternatively, if a ther- 
mal pulse occurs as the star evolves off the AGB, much of the 
remaining hydrogen is incinerated and nuclear processed ma- 
terial from deeper layers (e.g. C and O) is dredged up to the 
surface leading to a hydrogen deficient object (e.g. Iben et al. 
1983, Herwig et al.[T995). 

A number of empirical results are broadly supportive of this 
hypothesis. For example, the observed number ratios of hydro- 
gen rich to hydrogen-deficient central stars of planetary neb- 
ulae (2:1; Mendez fTWT) and DAs to DOs at T eff > 40000K 
(7:1; Fleming et al. 11986l > are roughly consistent with theoret- 
ical expectations (Iben & Tutukov [l984l l. The observed trans- 
formation of the hydrogen deficient PG1 159 central star of the 
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planetary nebula, Lo4, into a Wolf Rayet object and back over 
a period of months (Werner et al. 119921 and similarities in the 
measured abundances of the elements C and O in the atmo- 
spheres of Wolf R ayet (C/He~ 0.2 and O/He- 0.05, Koesterke 
& Hamann [19971 and PG1159 stars (C/H e~ 0.1 - 0.6 and 
O/He- 0.005 - 0.1, Dreizler & Heber 1998 ) confirm an evolu- 
tionary link between these two classes of object. Furthermore, 
the abundance patterns observed in the atmospheres of the DO 
white dwarfs (120000 >T eff > 45000K, C/He~ 0.001 - 0.01, 
0/He< 0.001; Dreizler 1998 1 suggest they are the descendents 
of the PG1159 stars (180000 >T eff > 65000K), gravitational 
settling having reduced the abundances in the former. 

Despite detecting ~ 1 10 hot DA white dwarfs (T e jf > 
25000K), the extreme-ultraviolet surveys of the ROSAT WFC 
and EUVE unearthed only one new DO white dwarf, RE 
J0503-285. An analysis of IUE and Voyager data and the 
optical identification spectrum revealed an effective tempera- 
ture and surface gravity of T e g a 70000K and log g =7.5- 
8.0 respectively (Barstow et al. 1994 ». These values are con- 
sistent with more recent determinations utilising more re- 
fined synthetic spectra and/or higher S/N optical data (T e ff — 
72660+g|K, log g = 7.50+° j3, Barsto w et a l. 12170171 T eff = 
70000K, log g=7.50, Dreizler & Werner [T996l . 

Somewhat counter-intuitively to it's detection as an EUV 
source, of the 14 DOs studied in detail by Dreizler & Werner 
( IT9961 and Dreizler ( fT^9"8l . RE J0503-318 is one of the most 
metal rich, with C/He* 0.005, O/He* 0.0005, N/He* 10 5 and 
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Table 1. Summary details of RE J0503-285, including coordinates and near-IR magnitudes obtained from the 2MASS 
All-Sky Point Source Catalogue. The exposure times used for the acquisition of the JH and HK near-IR spectra with the 
NTT and SOFI are also listed. 



Identity 


Name 


RA Dec J 


H 


Ks t exp 


(sees) 






12000.0 




IH 


HK 


WD0501-289 


RE J0503-285 


05 03 55.51 -28 54 34.6 14.65 +0.03 


14.77 + 0.06 


14.95 ±0.14 1080 


2880 



Si/Hew 10~ 5 . Nickel has also been detected in the photosphere 
at an abundance of Ni/H« 10~ 5 (Barstow et al. 2000 ). However, 
no convincing evidence of the presence of iron is found in any 
of the spectral datasets obtained to date, setting an upper limit 
on its abundance of Fe/H< 10~ 6 . To an extent this is an in- 
triguing result when one notes that the cosmic abundance ratio 
of Fe:Ni is 18:1, that Fe is observed to be more abundant than 
Ni in the atmospheres of hot DA white dwarfs by factors be- 
tween 1 to 20 and that self-consistent model calculations tak- 
ing into account gravitational settling and radiative levitation 
predict an apparent Fe:Ni ratio of ~1:3 (e.g. Dreizler 1998 1. 
Barstow & Sion ( 1994) report evidence of episodic massloss 
in a series of IUE spectra of RE J0503-285. It is plausible that 
the observed pattern of abundances may be reflective of this 
process. Alternatively, the abundance pattern may be related 
to the hypothetical dredging up of nuclear processed material 
from deeper layers caused by the late-thermal pulse, as has 
been proposed to explain the Fe deficiency observed in a num- 
ber of PG1159 stars (e.g. Miksa et al. 120021 . However, as the 
mass of RE J0503-285 maybe somewhat lower than the canon- 
ical DO value of 0.59 + 0.08M G (Dreizler & Werner [19961 it 
is a possibility that binary evolution (e.g. Vennes et al. I1998> 
has influenced or continues to interfere with the photospheric 
composition. Indeed, as discussed by Bragaglia et al. ( 1990) it 
is possible that close binary evolution produces preferentially 
non-DA white dwarfs with hydrogen-deficient envelopes. 

The lack of any prior published detailed examination of this 
particularly interesting object at near-IR wavelengths has mo- 
tivated us to obtain, during a recent study of DA white dwarfs, 
a JHK spectrum of RE J0503-285 to investigate external fac- 
tors possibly influencing its evolution and photospheric com- 
position e.g. the presence of a cool low mass companion. We 
present here the results of our analysis of this data. 

2. Observations 

2.1. Data acquisition 

A low resolution near-IR spectrum of RE J0503-285 was ob- 
tained using the ESO New Technology Telescope (NTT) and 
the Son-of-Isaac (SOFI) infrared instrument on 2003/12/10. 
The sky conditions at the La Silla site were good on this night 
with seeing typically in the range 0.6"- 1.0". We used the low 
resolution spectroscopic mode {A/6A ~ 950 with the 0.6" slit), 
in which coverage of the wavelength ranges 0.95 - 1 .64/um and 
1.53 - 2.52/mi is provided by the "blue" and the "red" grism 
respectively. The observations were undertaken using the stan- 
dard technique of nodding our point source target back and 
forth along the spectrograph slit in an ABBA pattern. The total 
integration times used for the blue and the red grism setup were 



1080s and 2880s respectively. To facilitate the removal of tel- 
luric features from the target spectra and to provide an approxi- 
mate flux calibration, a standard star (HIP28999) was observed 
immediately after the science integrations. This was carefully 
chosen to lie within ~ 0. 1 airmasses of the white dwarf. 

2.2. Data reduction 

We have used software routines in the STARLINK packages 
KAPPA and FIGARO to apply standard reduction techniques 
to our data. Further details of these procedures are given in 
Dobbie et al. (2005). Here we re-iterate that any features in- 
trinsic to the energy distribution of the standard star were iden- 
tified by reference to a near-IR spectral atlas of fundamental 
MK standards (Wallace et al. l2000l Meyer et al. [T998l Wallace 
& Hinkle 119971 and were removed by linearly interpolating 
over them. Furthermore, the flux levels of the data were scaled 
to (1) achieve the best possible agreement between the blue 
and the red spectrum of the white dwarf in the overlap region 
between 1.53 - 1.64/mi and (2) obtain the best possible agree- 
ment between the spectral data and the J, H and Ks photometric 
fluxes for RE J0503-285 derived from the 2MASS All Sky Data 
Release Point Source Catalogue magnitudes (Skrutskie et al. 
1997) where zero magnitude fluxes were taken from Zombeck 
( 1990). The reduced spectrum and 2MASS fluxes are shown in 
Figure 1. 

3. Analysis of the data 

3.1. Model DO white dwarf spectra 

We have generated a He+C+H synthetic white dwarf spec- 
trum with abundances C/He= 0.005 and H/He= 10" 5 , at the 
effective temperature and surface gravity shown in Table|2] We 
have used the latest versions of the plane-parallel, hydrostatic, 
non-local thermodynamic equilibrium (non-LTE) atmosphere 
and spectral synthesis codes TLUSTY (v200; Hubeny 1988 
Hubeny & Lanz [T995T ) and SYNSPEC (v48; Hubeny, I. and 
Lanz, T, private communication). The calculation included a 
full treatment of line blanketing and used state-of-the-art model 
atoms. In brief, the Hell ion incorporated the 19 lowest en- 
ergy levels, where the dissolution of the higher lying levels was 
treated by means of the occupation probability formalism of 
Hummer & Mihalas ( 1988), generalised to the non-LTE situ- 
ation by Hubeny, Hummer & Lanz ( 1994). Lines originating 
from the four lowest levels were treated by means of an ap- 
proximate Stark profile (Hubeny et al. 1994) during the cal- 
culation of the model structure and in the spectral synthesis 
step. The CIV ion incorporated levels up to n=14 where lev- 
els 9 < n < 14 were each represented by a superlevel. Where 
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Fig. 1. Near-IR spectroscopy (solid grey line) and 2MASS JHK photometry (filled circles) of the white dwarf RE J0503-285. A 
He+C+H non-LTE synthetic white dwarf spectrum of appropriate effective temperature, surface gravity and normalisation (solid 
black line) and a hybrid white dwarf + late-type dwarf model representing our estimated limit on the spectral type of a putative 
spatially unresolved companion (dotted black line, WD+M8) are overplotted. We have labelled the more prominent white dwarf 
Hell lines and telluric water vapour features present in this dataset. 



Table 2. Summary of additional physical parameters of RE 
J0503-285 taken from Barstow et al. (1994, 2000). Shown 
in brackets are the values of effective temperature and sur- 
face gravity, lying within the uncertainties of Barstow et al's 
analysis, which place RE J0503-285 at the largest likely dis- 
tance. We adopt this latter value in setting a limit on the spec- 
tral type of a putative cool companion. Note, the distances 
have been estimated using pure-C core evolutionary models 
of Wood ( 1995 ) which exclude a hydrogen layer. 



Teff(K) 


logg V 


D(pc) 


72660(75613) 


7.50(7.35) 13.90 ± 0.01 


152(176) 



available, photoionization cross sections were obtained from 
TOPbase. Transitions between the lowest level and those with 
n< 5 were represented by Stark profiles. 

The synthetic energy distribution has been normalised to 
the V magnitude of RE J0503-285 and convolved with a 
Gaussian to match the resolution of the SOFI spectrum. This 
is shown overplotted on the observed data in Figure 1. It is 
worth noting here that our modelling indicates that at the effec- 
tive temperature of this white dwarf, the colours V-K, J-H and 
H-K are rather weak functions of T e u . 

3.2. The search for a cool companion 

We have examined Figure 1 for significant differences in the 
overall shape or level between the observed and synthetic 
fluxes which can be consistent with the presence of a cool com- 
panion. Further, we have searched for specific features in the 
spectrum typical of the energy distributions of M or L dwarfs 
e.g. K I and Na I absorption at 1.25/mi and 2.20yum respec- 



tively, CO at 2.3/mi and H2O centred on 1.15, 1.4 and 1.9/mi. 
However, no convincing evidence for such has been found. 
Hence, we have added empirical models for low mass stellar 
and substellar objects, full details of which are given in Dobbie 
et al. (2005 1, to the white dwarf synthetic spectrum and com- 
pared these composites to the IR data to set an approximate 
limit on the spectral type of a putative cool companion. 

The fluxes of the empirical models have been scaled to a 
level appropriate to a location at d=10pc using the 2MASS J 
magnitude of each late-type object and the polynomial fits of 
Dahn et al. (2002 1 and Tinney et al. d2003l to the Mj versus 
spectral type for M6-M9 and L0-L8 field dwarfs/brown dwarfs 
respectively. These fluxes have been further reduced by a factor 
1 .4, corresponding to the rms dispersion in the Mj versus spec- 
tral type relationship of Tinney et al. (2003 1. Subsequently, the 
fluxes have been re-calibrated to be consistent with the distance 
of the DO as derived from the measured V magnitude, the ef- 
fective temperature and theoretical My and radius from evolu- 
tionary models of pure-C core white dwarfs which include only 
a He layer mass of 1O~ 4 M (Wood 1995 1. As the effective tem- 
perature and surface gravity and hence our distance determi- 
nation of RE J0503-285 are considerably less well constrained 
than those of the DA white dwarfs we have previously investi- 
gated (e.g. Dobbie et al. 2005J, for a conservative limit we have 
adopted the largest distance estimate consistent with Barstow 
et al.'s analysis and the error limits given therein (see Table 2). 
In fact the neglect of line blanketing arising from metals other 
than carbon present in the atmosphere of REJ0503-285 impacts 
on the temperature estimate in a such a way as to reinforce 
our cautious approach (e.g. see Barstow et al. 2000 1. Starting 
with L8 we have progressively added earlier spectral types to 
the synthetic white dwarf spectrum, until it could be concluded 
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with reasonable certainty that the presence of a companion of 
that effective temperature or greater would have been obvious 
from our data, given the S/N. 



4. Results 

4.1. A limit on the mass of a pututive companion to RE 
J0503-285 

The effective temperatures of very-low-mass stars and substel- 
lar objects remain a function of both mass and age on the 
timescale typical of the lifetime of the white dwarfs likely pro- 
genitor (~Gyrs). Although we can estimate the cooling age of 
RE J0503-285 using theoretical evolutionary models, as we 
don't know with any certainty the mass and hence lifetime of 
its progenitor star (the initial-mass final-mass relationship is 
poorly constrained for white dwarfs of this mass and in any 
case may not apply to DOs) a robust estimate of the age of a 
putative associate of the DO is not possible. Therefore, to use 
our limit on spectral-type to constrain mass, we instead assume 
a range of ages broadly encompassing the likely value. We as- 
sign an approximate effective temperature to the spectral type 
limit shown in Table [3] using the polynomial fit described in 
Table 4 of Golimowski et al. (2004). Subsequently, we refer 
to the low mass stellar/substellar evolutionary models for solar 
metallicity of Baraffe et al. ( 2003 1, using cubic splines to inter- 
polate between their points, to estimate corresponding mass at 
ages lGyr, 5Gyrs and lOGyrs as shown in Table 3. 

We note that there are a number of examples of post com- 
mon envelope white dwarf + main sequence binaries, where a 
low mass companion has survived a phase of common envelope 
evolution e.g. RE J0720-318, REJ1016-05, and RE J2013+400 
and now interacts with the white dwarf to influence the com- 
position of its photosphere. Our result argues strongly against 
the presence of an unresolved cool companion to this white 
dwarf with M> O.O85M and we can now rule out that the pat- 
tern of abundances observed in the atmosphere of RE J0503- 
285 is influenced by ongoing interaction with such an object. 
However, our result does not exclude the possibility that the 
evolution of RE J0503-285 was affected by the interaction 
of the white dwarf's progenitor star and a close companion 
(M> 0.085M o ) during a common envelope phase from which 
the latter failed to emerge. Further, there remains the possibility 
of a close companion with a mass below our detection thresh- 
old. This may seem unlikely given that an extensive near-IR 
study of a sample of 371 white dwarfs, summarised by Farihi, 
Becklin & Zuckerman (2005), has revealed only two secon- 
daries with spectral types later than M, GD165B (L4; Becklin 
& Zuckerman 1998) and GD1400B (L6; Farihi & Christopher 
2004, Dobbie et al. 2005), the former at a separation a~120AU 
from its white dwarf primary. Nevertheless, there are persistent 
although perhaps controversial claims for L dwarf secondaries 
in some cataclysmic variables such as EF Eri (e.g. Howell & 
Ciardi 2001, Harrison et al. 2003). 



Table 3. Limiting spectral type, temperature and mass of 
a putative cool companion to RE J0503-285. The effective 
temperature is estimated from the polynomial relation de- 
tailed in Table 4 of Golimowski et al. (2004 1. Furthermore, 
we provide rough upper limits on the mass of a putative cool 
companion as a function of age, by comparing this effective 
temperature to the predictions of the evolutionary models of 
Baraffe et al. d2003l . 



SpT T cff (K) 


lGyr(M G ) 


5Gyr(M G ) 


lOGyr(M ) 


M8 2500 


0.084 


0.085 


0.085 



4.2. Possible future work 

We point out that RE J0503-285 and the other DO white dwarfs 
are potential targets for Spitzer imaging to search for a mid- 
IR excess due to companions cooler than spectral type M. Our 
non-LTE model of REJ 0503-285 indicates the white dwarf flux 
at 10;um to be ^0.04mJy. From the Spitzer spectrum of DENIS- 
P J0255-4700 (Roellig et al. 2004J, which we determine to re- 
side at d=5.3pc using the 2MASS magnitude (J=13.25+0.03) 
and the Mj versus spectral type relationship of Tinney et al. 
(2003), allowing for a 25% uncertainty in the absolute calibra- 
tion of the mid-IR data and the uncertainties in late-type dwarf 
fluxes and the white dwarf distance discussed in Section 3.2, we 
estimate the flux level of an L8 companion to RE J0503-285 to 
be x;0.004mJy. As the absolute accuracy of Spitzer photome- 
try is expected to be ~5%, it should be possible to extend this 
search to the regime of the coolest L dwarfs, corresponding to 
substellar masses. Further, improved line broadening theories 
specifically tailored to helium plasmas would offer the possi- 
bility of reducing the uncertainties on our estimates of the ef- 
fective temperatures and surface gravities of DO white dwarfs 
and allow more robust determinations of their distances. 
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